Direct observation of current in type I ELM filaments on ASDEX Upgrade 
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Magnetically confined plasmas are often subject to relaxation oscillations accompanied by large transport 
events. This is particularly the case for the high confinement regime of tokamaks where these events are termed 
edge localized modes (ELMs). They result in the temporary breakdown of the high confinement and lead to high 
power loads on plasma facing components. Present theories of ELM generation rely on a combined effect of 
edge current and the edge pressure gradients which result in intermediate mode number (n = 10 — 15) structures 
(filaments) localized in the perpendicular plane and extended along the field line. It is shown here by detailed 
localized measurements of the magnetic field perturbation associated to an individual type I ELM filament that 
these filaments carry a substantial current. 
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PACS numbers: 52.35.PY 52.55.Fa 52.55.RK 52.70.Ds 

Edge localised modes (ELMs) are short (ms) breakdowns 
of the high confinement regime (H-mode), which is envis- 
aged to be used in power producing tokamaks. Due to the 
high power fluxes associated, the presence of ELMs poses de- 
mands on the design of plasma facing components, that are 
hard to meet. Thus understanding, with the aim of control- 
ling, ELM events is one of the foremost priorities in fusion 
research. Moreover, apart from the interest for fusion oriented 
plasmas, the interest in ELM physics is enhanced by some fas- 
cinating analogies with sporadic explosive events as observed 
for example in solar flares fl]] or in magnetic substorms yfl. 

ELMs are at present thought to originate from a combina- 
tion of current and pressure gradient driven MHD modes J^, 
and are known to result in a medium number (n w 10 — 15) 
of structures 0-01, localised in the perpendicular plane and 
extended along the magnetic field. These filaments travel 
through the Scrape Off Layer (SOL) and have been measured 
using Langmuir probes on various machines, see e.g. |8y, and 
observed through high speed cameras |01 or Gas Puff Imaging 
diagnostic jjj. In this letter we present results of direct mea- 
surements of all three components of the magnetic field per- 
turbations associated to ELM filament structures in the SOL 
together with an estimate of the current carried by filaments. 

Magnetic fluctuations associated with ELMs are usually 
believed to originate from MHD activity. Measurements of 
the magnetic activity have previously been performed using 
magnetic pickup coils close to the vessel wall J3, fioll or on 
insertable probes 111 ill . These measurements generally take 
place far from the filaments in comparison to their radial ex- 
tent. This makes it difficult to observe the magnetic pertur- 
bation going along with individual filaments and to examine 
the magnetic fine structure of the ELMs. Both would result 
in important information such as the excursion of magnetic 
field lines from their equilibrium position and if the ELMs are 
associated with reconnection events (lj] . 

Here we report on measurements of type I ELMs carried 
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FIG. 1. Top: Degree of Polarization (DOP) analysis as a function 
of time and frequency. Bottom: Ion saturation current (I 5o; ). Right: 
Coherence between \ ml and b p computed in the ELM and inter-ELM 
phases. The time intervals are shown in the bottom panel with col- 
ored boxes. 



out at ASDEX Upgrade (AUG) tokamak by means of a newly 
constructed probe head described in detail elsewhere IU2I1 . 
The probe head consists of a cylindrical graphite case (jzf=60 
mm) which holds six graphite pins. One of the tips measured 
the ion saturation current, one tip was swept, whereas all the 
others had been kept floating. Inside the case, 20 mm behind 
the front side, a magnetic sensor measuring the time deriva- 
tive of the three components of the magnetic field is mounted. 
The sensor has a measured bandwidth of 1 MHz. A similar 
probe with combined electrostatic and magnetic signals has 
already been used in AUG II ill , but with a wider magnetic 
coil and a larger radial separation between electrostatic and 
magnetic sensors. The probe head, mounted on the fast recip- 
rocating midplane manipulator, is inserted from the low field 
side of the torus for 100 ms 12 mm inside the limiter posi- 
tion. The data have been obtained in type-I ELMy plasma dis- 
charges (# 23158,23159,23160,23161,23163) with a toroidal 
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FIG. 2. (a)Ion saturation current (b) Poloidal and radial component 
of the magnetic field. The red box highlights the time interval when 
the two components change their phase relation, (c) Hodogram of 
the radial/poloidal magnetic field component. The closed loop in red 
refers to the time intervals highlighted in panel (a,b). 
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magnetic field of -2.5 T, 0.8 MA of plasma current, 6.5 x 10 
m~ 3 central electron density and a ^95 value of 5.2. 

During the analysis we have used the ion saturation cur- 
rent to infer the passing of a type-I ELM filament structure in 
front of the probe, as done for example in 1 1311 . In analyzing 
the ELM events we employ the idea that the magnetic sig- 
nal during the ELM can be separated into different frequency 
components. Higher frequencies, above few hundred kHz, we 
presume to be generated mostly by Alfvenic activity or high 
frequency turbulence. These are not considered in the present 
analysis, also because of the frequency cut-off by the graphite 
shielding 1 1 211 - Additional MHD activity will still be present 
in the signal at frequencies below 20 kHz, but during ELM fil- 
aments most of the signal is presumed to originate from slowly 
varying currents convected with the filaments. This is justified 
by the so-called Degree of Polarization (DOP) analysis 11411 . 
which is a test for a plane wave ansatz, quantifying how well 
the relation k ■ B = is satisfied. It is based on the evalu- 
ation and diagonalization of the spectral matrix S = (B*Bj), 
calculated in Fourier space. The DOP represents a measure to 
determine whether S represents a pure state quantifying how 
much one single eigenvector approximates the state. A high 
value of DOP implies that the fluctuations considered are co- 
herent over several wavelengths, and thus the method can be 
used to distinguish between propagating modes and coherent 
localized fluctuations. Figure [U panel (a), shows the results 
of DOP analysis as a function of frequency and time. The 
temporal evolution of the ion saturation current is depicted 
in the lower panel. A sudden drop in the DOP is observed 
in correspondance with a steep increase of I sat signals: this 
implies that the magnetic field fluctuations during an ELM 
can be better represented by coherent structures than by plane 
wavepackets. A coherence analysis between the ion satura- 
tion current and the poloidal component of the magnetic field 
(Fig. [U(c)) reveals an increase of coherence between the two 
signals during the ELM activity. A closer look on the be- 
havior of magnetic fluctuations associated to an ELM event is 
shown in Figure [2] (b). When the magnetic fluctuation ampli- 
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FIG. 3. (a)Trajectory of ELM filament associated magnetic field ex- 
cursions in all three spatial directions. The direction of minimum 
variance is shown (black line) together with the direction of the equi- 
librium field (blue line), (b) Hodogram of the magnetic perturbation 
associated with ELM current filament reconstructed in the maximum 
variance plane. The red line shows the elliptical fit. 



tude increases, the radial and poloidal components (b r and b p 
) change their phase relation as highlighted by the color box. 
We now consider the hodogram of the perpendicular magnetic 
perturbation (Fig. 0(c)), i.e. the magnetic field perturbation 
trajectory in the b r — b p plane, during the passing of the ELM. 
Here one closed orbit corresponds to the time interval marked 
in Figure [2] b. Closed loops are compatible with the passing 
of current filaments in front of the probe. Outside the ELM, 
the magnetic field perturbation exhibits an almost linear po- 
larization in the perpendicular plane. This shows clearly that 
magnetic activity in between ELMs (wavelike) differs qualita- 
tively from the magnetic perturbation during an ELM, which 
seems to be due to the motion of current filaments. 

The measurements of all three components of the magnetic 
perturbation allows to check the alignment of the current fila- 
ments directly. In Figure [3] (a) the trajectory of the magnetic 
field excursion during an ELM event in all three magnetic field 
components is shown. The time interval selected corresponds 
to the one highlighted in Figure|2] A closed elliptical loop ly- 
ing in a plane slightly tilted with respect to the local frame of 
reference is observed, as expected for filamentary structures. 
The direction normal to this plane can be determined by us- 
ing the minimum variance analysis 1 15]. The method is based 



on the solution of the eigenvalue problem £v=i M" v ii v 
where M^ v = (B p B v ) — (B p )(B v ) is the Magnetic Variance 
Matrix, the brackets indicating the mean values averaged over 
the time the structure spends traveling in front of the probe, 
/j,V = 1,2,3 denoting the cartesian components of theX,Y and 
Z of the chosen system, A, and n being respectively the eigen- 
values and eigenvectors of the system. The eigenvectors cor- 
responding to the eigenvalues X\, X2 and A3, represent the di- 
rection of maximum, intermediate and minimum variance of 
the magnetic field respectively. The eigenvector correspond- 
ing to the minimum eigenvalue corresponds to the direction 
of minimum variance and is normal to the plane spanned by 
the magnetic field perturbations. This direction is observed 
to be parallel to that of the equilibrium magnetic field, shown 
in the same plot by a blue line. This confirms the hypothesis 
that the ELM filament is aligned with the equilibrium mag- 
netic field. From the sense of the polarization the direction 
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FIG. 4. Hodograms calculated for a monopolar current distribution 
(9W ) and bipolar one ($). 



of the current is found to be colinear with the plasma current. 
The method allows the determination of the rotation matrix, 
so that the hodogram in the plane perpendicular to the current 
filaments may be reconstructed. This hodogram is shown in 
Figure [3] (b) together with a fit to an ellipse. The shape of 
the hodogram calculated in the rotated frame of reference can 
be used to determine the type of current distribution associ- 
ated to an ELM. Theories have b oth p roposed ELM filaments 
to be associated to monopolar Ji, 16] and bipolar current dis- 
tributions Il7tl . Up to now no clear evidence of one mecha- 
nism predominant with respect to the other has been reported. 
Figure |4] shows the anticipated shape of the hodogram in the 
bipolar and monopolar cases. The hodogram in the bipolar 
case exhibits a cardioid-like shape with a cusp at the origin, 
which represents a distinct signature and which can not be 
recognized in the experimental data (see Fig [3] (b)). To re- 
inforce the statement we note that a bipolar-like hodogram 
with the presence of a cusp has been previously recognized 
in lfl8ll (cfr. fig 3(b)) where it has been associated to a direct 
measurement of a bipolar current structures. We emphasize, 
however, that the filaments observed in ll 811 are of a com- 
pletely different origin as induced by Drift- Alfven turbulence 
11911 . A further corroboration of the hypothesis of monopo- 
lar current filaments results from considering the quadrants 
occupied by the magnetic field trajectory. Indeed while the 
monopolar current hodogram regularly occupies two of the 
quadrants, the bipolar one always spans three of the quad- 
rants: this latter behavior is not observed in the experimental 
data as can be verified by comparing Figure [3] (b) and Figure 
[4] These results make us confident that the magnetic fluctu- 
ations observed are indeed generated by a monopolar current 
distribution. Under this basic assumption the current carried 
by this filament may be estimated. In the rotated frame of ref- 
erence for a circular monopolar symmetric filament drifting 
in the X\ direction, the two perpendicular magnetic field com- 
ponents may be written as bi, = ~'J> B °" and fcv, = 

where Bq = and a is the distance between the trajectory 



of the center of the filaments and the probe, representing also 
the distance of closest approach, ro is the radius of the fila- 
ments and Iq its current. The distance a can be approximated, 
assuming the filament to propagate with a constant velocity 
along the X\ direction: a = Atvxj where At is the time delay 
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FIG. 5. (a) Time traces of the two perpendicular components of mag- 
netic field during an ELM filament, (b) Total perpendicular magnetic 
field as a function of time. Superimposed the results of a fit for an 
expected monopolar (red) and bipolar (blue) current distribution both 
modulated by an exponential decay. 



between the maximum of b%^ and the maximum/minimum of 
bx 2 , where bx l and bx 9 are the projections of the magnetic field 
in the maximum-intermediate variance plane. Thus within 
these approximations the current may be estimated noting that 



_ l mk 

2 2xa 



4xAtv x 



Applied to the experimental 



data, this computation is equivalent to calculate the quantities 
depicted in Figure [5] where At has been calculated both at the 
maximum and at the minimum of bx r The two values At\ 
and A/2 are equal to 38 /js and 42 /js respectively. The esti- 
mate of the current relies on the knowledge of local velocity 
in the X\ direction. T his p ropagation has been reported for 
ASDEX Upgrade, e.g. II ll 12011 . The experimental setup does 
not allow a reliable local measurements of v, or v p , neither 
can we rely in the presented shots on correlation analysis of 
wall mounted probes (which are 160 degrees separated from 
the manipulator in the toroidal direction). Most likely these 
signals are indeed dominated by MHD mode activity in the 
confined plasma rather than current filaments in the SOL. As 
the best approximation we thus assume as radial propagation 
the most probable value v r = 1.2 km/s as determined in lEoll 
and calculate \>x 1 = v,-/ cos(Z(A,i — r)), where Z(k\ — r) is the 
nagle between X\ and the radial direction which for the present 
case is of the order of 7°. From If20ll we also estimate a stan- 
dard deviation of v r as o Vr = 700 m/s, which ensure that 7 1 % 
of the events are observed within (v r ±a Vr ). With these val- 
ues the average distance of closest approach is approximately 
4 cm and using the average value between the minimum and 
maximum of bx~ we obtain an estimate of 1.9 kA. The to- 
tal perpendicular field yjb^ + b^\ 2 is shown in figure [5] (b). 
The magnetic field variation can be fitted with a function of 



the form / = 



-('-fo)/T< 



, to being the time 



instant corresponding to the maximum of the total perpendic- 
ular field and i e the e-folding time, and assuming the mag- 
netic field is the response from a passing monopolar current 
filament. The good quality of the fit provides additional sup- 
port for the monopolar nature of the filament, also comparing 
with the fit expected from a bipolar current distribution shown 
in blue line which exhibits an higher % 2 . The average decay 
time is determined as % e w 200 /js. In order to increase the 
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FIG. 6. (a) Joint PDF of At and (X\ = a) (b) )Histogram of the 
angle between \\ and the radial direction, (c) Filaments current es- 
timate Vi distance from the separatrix estimated with v r = 1.2 km/s. 
Error bars results from the propagation of velocity uncertainty in the 
estimate of a and /. 



statistic reliability of the previous estimate we have analyzed 
events from 5 different shots with similar conditions. The re- 
sults are shown in Figure [6] In the panel (a) the joint PDF 
of the independent experimental values At and bx 2 used for 
the current evaluation is shown, highlighting how the bulk of 
the filaments have values of approximately 25 /js and 5 mT. 
In the panel (b) the histogram of the angle between Xi and 
the radial direction is shown, showing that the A-i direction 
is a complex combination of radial and poloidal propagation. 
Finally in the panel (c) the current estimated accordingly to 
the aforementioned formula is plotted vs the position of the 
filaments with respect to the separatrix, taking into account 
the position of the manipulator and the distance of closest ap- 
proach estimated. The errors shown represent the influence 
of the velocity uncertainty on the estimate of a and conse- 
quently on /. The bulk of the distribution of these filaments 
are thus observed in the SOL, even taking into account the 
uncertainty on v r : the median of the distribution of the fila- 
ments detected in the SOL within the error in their position is 
1.4 kA, corresponding to a y'|| rj 4.5MA/m 2 for 1 cm radius 
filaments. These values are consistent with measurements of 
edge current as seen for example in 12 111 . It has already been 



supposed in 12211 that, in case ELM breakdown is driven by 
peeling instabilities, this will lead to a flattening of the cur- 
rent density profile over the separatrix and a subsequent in- 
crease of currents flowing into the SOL. These currents will 
be nearly force-free and will be accompanied by poloidal halo 
currents closing through the divertor tiles. Actually these cur- 
rents have been measured 112l 12311 with values up to few tens 
of kA, supporting the presence of rapid flow of toroidal cur- 
rents from the plasma into the SOL The resulting histogram 
may also be compared with the estimate for ELM filaments 
in JET given in [24]: in this case the most probable current 
has been found of the order of 450 A, but with a lower ra- 
dial velocity postulated. For completeness it must be noted 



that the same current density was estimated in 11 111 assuming 
the magnetic perturbation to be induced by a rotating helical 
structure with a bi-directional current close to, but still inside, 
the separatrix. The value found is higher than the measured 
jsat current density to the ion-biased tip (approximately 50 
kA/m 2 ). Concluding in this letter we have provided evidence 
that ELM filaments carry considerable currents for which we 
have found a reasonable estimate. The magnetic signals dur- 
ing ELM filaments differ substantially from wave activity in 
between ELMs. We have shown that the current in the fila- 
ments is co-aligned to the plasma current and approximately 
of a magnitude as expected for the edge. The current flows 
along the unperturbed magnetic field lines and has a unidi- 
rectional nature. This poses the question where the filament 
currents close in the SOL and why such high current densi- 
ties are sustained in the ELM filaments. We hope that future 
experiments will contribute to answer these questions, which 
will throw new light on the instability mechanisms for ELMs. 
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